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CALCULATING  BOREHOLE  GEOMETRY  FROM  STANDARD  MEASUREMENTS 
OF  BOREHOLE  INCLINOMETRY 

Kenneth  C.  Jezek  and  Richard  B.  Alley 

INTRODUCTION 

The  most  common  method  for  estimating  borehole  geometry  in  glaciers  is 
standard  borehole  inclinometry.  Inclinometry  data  collected  during  the 
drilling  operation  are  used  to  monitor  the  progress  of  the  drill  (Gundestrup 
et  al.  1982)  and  eventually  are  used  to  orient  and  interpret  the  stratigraphy 
in  the  core.  Subsequent  remeasurements  of  borehole  geometry  are  the  most  ac¬ 
curate  methods  for  estimating  intraglacial  stjrain  rates  (Savage  and  Patterson 
1963,  Garfield  and  Ueda  1976). 

Inclinometry  logs  measure  local  changes  in  the  inclination  and  azimuth 
of  the  borehole.  A  method  successfully  used  by  Gundestrup  and  Hansen  (in 
press)  to  log  the  deep-drill  hole  at  DYE-3  on  the  Greenland  Ice  Sheet  employs 
two  pendulum  type  inclinometers  with  their  sensitive  axes  orthogonal.  A  mag¬ 
netic  compass  is  used  to  determine  the  magnetic  azimuth  of  the  sensitive 
axes.  Measurements  are  made  at  discrete  points  spaced  on  the  order  of  tens 
of  meters  along  the  borehole  and  are  converted  to  values  of  borehole  inclina¬ 
tion,  azimuth  and  arc  length.  (Nearly  continuous  yet  rapid  logging  is  possi¬ 
ble  with  a  new  technique  that  tracks  a  radar  target  lowered  into  the  borehole 
[Jezek,  in  press],  but  the  method  is  not  yet  as  accurate  as  the  standard  pro¬ 
cedure.)  The  precision  and  reproducibility  of  these  measurements  are  very 
good  —  1%  of  the  measured  inclination  and  -  2°  in  the  azimuth*  —  but  a 
dilemma  arises  in  choosing  an  analytical  method  for  computing  the  coordinates 
by  integrating  the  inclination  and  azimuth  data.  Several  methods  can  be 
found  in  the  literature  of  the  petroleum  industry  (Rivero  1971,  Blythe  1973, 
Taylor  and  Mason  1971). 

In  this  paper  we  will  describe  the  details  of  a  method  we  used  in  ana¬ 
lyzing  inclinometry  data  collected  in  the  borehole  at  DYE-3,  Greenland.  Our 
model  is  similar  to  that  proposed  by  Walstrom  et  al.  (1972)  (their  model  4) 

*Personal  communication  with  B.L.  Hansen,  of  the  Polar  Ice  Coring  Office, 
Lincoln,  Nebraska,  1984. 


and  makes  Che  simple  assumption  that  inclination  and  azimuth  vary  linearly 
with  arc  length  between  measuring  points.  Our  model  suffers  from  a  universal 
flaw  in  these  approaches,  which  is  the  ad  hoc  assumption  about  how  the  bore¬ 
hole  varies  between  measuring  points.  However,  we  claim  that  our  method  can 
be  used  to  create  physically  reasonable  models  even  in  cases  where  the  bore¬ 
hole  is  believed  to  behave  in  a  complex  fashion.  Furthermore,  we  feel  our 
technique  is  convenient  and  straightforward  to  understand. 

METHOD 

We  assume  that  the  inclination  0  and  azimuth  a  are  different  linear 
functions  of  arc  length  S  between  measurement  points. 


G  =  aS  + 

©0 

(la) 

a  =  cS  + 

«o 

(lb) 

where 
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and  where  0o  and  aQ  are  the  inclination  and  azimuth  at  the  top  of  an  interval 
of  arc  length  As,  and  ©i  and  04  are  the  inclination  and  azimuth  at  the  bottom 
of  that  interval.  We  imagine  the  borehole  as  a  smooth,  three-dimensional 
curve  and  describe  it  by  a  set  of  three  equations  parameterized  against  arc 
length  (Fig.  1) , 
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Now  it  is  straightforward  to  substitute  eq  1  into  eq  3  and  integrate. 
The  result  is  that  the  Cartesian  coordinates  x,  y  and  z  at  the  bottom  of  the 
interval  of  arc  length  are  given  by 

x  =  x'  -  P( A  +  B)  +  Q(E  +  F)  -  R(E  -  F)  -  T(A  -B) 

+  p  ( - ® j  f - £ - ] 

r  (a-c)  (a+c)  J  1  (a-c)(a+c) 


(4a) 


Figure  1.  Geometry  for  computing  the 
coordinate  transformation.  The  x-y 
plane  at  z  =  0  represents  the  surface 
of  the  glacier.  The  heavy,  curved 
line  represents  the  borehole.  The 
azimuthal  angle  a  is  measured  relative 
to  the  positive  x  direction  (which  is 
true  north  for  the  coordinates  of  the 
DYE-3  borehole  data).  The  inclination 
angle  6  is  measured  relative  to  the 
z-axis. 
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and  where  x' ,  y'  and  z'  are  the  Cartesian  coordinates  at  the  top  of  the  in¬ 
terval.  The  equations  for  x  and  y  are  different  for  the  case  where  a  =  c; 
then 


x  =  x*  +  -1--2-^-as-)  (P  -  T)  +  ~  (Q  "  R>  + 
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The  equation  for  z  is  different  for  the  case  when  ©Q  =  Oi  and  a  =  0;  then 


z  *  z*  +  S  cos  0o. 


(8) 


This  model,  like  most  other  published  models,  assumes  that  all  informa¬ 
tion  about  borehole  orientation  along  an  interval  is  contained  in  data  meas¬ 
ured  at  the  endpoints  of  that  interval.  In  some  cases,  borehole  orientation 
may  be  related  to  the  orientation  in  the  previous  interval.  This  is  likely 
to  be  true  for  a  new  borehole  in  a  glacier  because  drilling  tends  to  progress 
along  a  smooth  curve,  but  may  be  less  likely  for  an  old  borehole  where  any 
shear  deformations  (which  could  be  caused  by  material  discontinuities  such  as 
dust)  will  distort  the  borehole  in  an  unpredictable  way.  If  the  deformation 
rate  is  small,  then  a  preferred  model  might  be  quadratic  expressions  for  0 
and  a  such  as 

0  =  eS2  +  fS  +  0O  (9  a) 

a  -  gS2  +  hS  +  aQ  (9  b) 

where 
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and  where  0_i  and  a_i  are  the  inclination  and  azimuth  at  the  point  above  the 
interval  in  question  and  where  AS  is  assumed  constant  for  the  length  of  the 
borehole.  Equation  9  can  be  substituted  into  eq  3  and  integrated  numerically 
to  trace  the  borehole.  (We  note  that  if  borehole  orientation  in  an  interval 
is  related  to  even  longer  lengths  of  borehole,  then  recursion  formulae  for 
higher  order  polynomials  can  be  generated  for  computing  the  transformation.) 

DISCUSSION  AND  CONCLUSIONS 

We  entered  the  equations  describing  linear  variations  of  0  and  a  on  a 
programmable  hand  calculator  and  used  a  simple  difference  test  to  choose 
whether  eq  4,  7  or  8  would  be  used  for  each  interval  (Appendix  A).  As  an 
example,  we  applied  the  method  to  the  logging  data  (Gundestrup  and  Hansen,  in 
press)  from  the  DYE-3  borehole  (Fig.  2).  We  averaged  data  where  multiple 
measurements  were  available  and  used  smooth  curves  to  extrapolate  0  and  a 
from  the  shallowest  measurements  at  the  100-m  depth  to  the  surface.  The  re¬ 
sults  are  shown  in  Figure  3. 

Wo  also  compared  the  linear-variation  and  quadratic-variation  models  for 
part  of  the  borehole  (Appendix  B).  The  Cartesian  coordinates  generated  by 
quadratic  variations  are  only  about  0.5%  different  from  those  generated  by 
the  linear-variation  model  (Table  1). 

A  drawback  to  the  linear-variation  model,  which  is  minimized  somewhat  in 
the  quadratic  model,  is  that  the  variation  of  0  and  a  along  the  entire  length 
of  the  borehole  is  continuous  but  not  smooth.  This,  however,  is  a  consider¬ 
able  improvement  over  other  methods  (such  as  assuming  that  the  angles  of  in¬ 
clination  and  azimuth  measured  at  the  bottom  of  the  interval  are  valid  over 


Table  1.  Comparison  of  linear-variation  and  three-point  quadratic-variation 
models  for  1900  to  1925-m  depth.  Run  with  n=6  in  quadratic  model. 


Difference 


s(m) 


Inclination 

6° 


Azimuth 

a° 


Ay(m) 


Linear 


_  r  Quad-Lin  >  x  100% 

Quadratic  Quad 


Figure  3.  Spatial  coordinates  of  the  DYE-3 
borehole  computed  by  applying  the  linear- 
variation  model  to  the  inclinometry  data  col¬ 
lected  by  Gundestrup  and  Hansen  (in  press). 

The  open  ellipses  represent  the  projection  of 
the  borehole  onto  the  surface  of  the  glacier. 
The  solid  circles  represent  points  along  the 
borehole  at  different  depths. 

the  entire  interval  or  using  average  values  for  the  up-segment  and  down- 
segment  coordinates  based  on  the  angles  measured  at  the  end  points  of  the 
segment)  which  produce  models  of  0  and  a  that  are  not  continuous. 

Our  purpose  has  been  to  describe  the  transformation  of  azimuth,  inclina¬ 
tion  and  arc-length  data  to  Cartesian  coordinates,  rather  than  to  examine  a 
particular  logging  experiment.  Thus,  we  limit  our  conclusions  to  the  state¬ 
ment  that  the  linear-variation  method  is  conceptually  simple,  straightforward 
to  implement,  and  that  it  _an  be  expanded  to  nonlinear  variations. 
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APPENDIX  A:  TEXAS  INSTRUMENTS  Ti-59  CALCULATOR 
PROGRAM  FOR  LINEAR-VARIATION  MODEL 

Memories  are  numbered  01  to  39. 

Keystrokes  preceded  by  the  shift  (2nd)  key  are  underlined. 

Abbreviations  used: 

INV  -  inverse  (of  the  following  keystroke). 

x  ^  t  -  puts  displayed  number  into  test  register. 

x>t  -  transfers  control  to  the  following  location  if  the  dis- 

played  value  (x)  is  greater  than  or  equal  to  the  value  i 
the  test  register  (t);  otherwise,  control  does  not  pass 
to  the  indicated  location  but  goes  to  the  next  statement 

2 

INV  x>t  Ex  -  will  pass  control  to  statements  following  label  E  if  x  i 

less  than  t;  otherwise,  it  will  square  x. 


cos 

— 

cosine. 

sin 

- 

sine. 

STO 

XX 

- 

stores  the  displayed  value  in  memory  xx. 

SUM 

XX 

- 

adds  the  displayed  value  to  the  value  in  memory  xx,  and 
stores  the  sum  in  memory  xx. 

RCL 

XX 

- 

recalls  the  value  in  memory  xx. 

GTO 

label 

- 

GO  TO  label,  unconditionally  transfers  control  to  the 
indicated  location. 

LBL 

- 

Label. 

OP  xx 

special  calculator  operation  (such  as  partitioning  the 
calculator  between  addressable  memories  and  program 
steps) . 

R/S 

- 

run/stop;  stops  operation  if  calculator  is  running. 

*  -  multiplication. 

4  0£  17  -  partitions  calculator  for  39  memories,  639  steps, 

puts  calculator  into  radians  mode. 


Rad 


Enter  program  into  calculator. 

Store  0O  (in  radians  =  degrees  *  it  t  180)  in  memory  02. 
Store  a0  (in  radians  =  degrees  *  u  *  180)  in  memory  04. 
Store  9i  (in  degrees)  in  memory  26. 

Store  ai  (in  degrees)  in  memory  27. 

Store  As  (in  meters)  in  memory  05. 


Store  0.0001  in  memory  30  —  used  in  test  to  decide  which  sets  of  equations  to 


Press  A. 

The  calculator  displays  Ax  (in  meters). 

The  calculator  contains  Ax,  Ay  and  Az  in  memories  20,  21  and  22  respectively. 

For  subsequent  intervals  -  enter  new  9 i  and  into  memories  26  and  27  (in 
degrees) : 

-  Enter  new  As  (this  is  only  necessary  if  As  changes) 

-  Press  A. 

After  each  calculation,  the  calculator  displays  Ax  for  that  interval: 

-  Memories  20,  21  and  22  contain  Ax,  Ay  and  Az ,  respectively,  for 
that  interval. 

-  Memories  23,  24  and  25  contain  x,  y  and  z  (EAx,  EAy  and  EAz) , 
respectively,  to  the  bottom  of  that  interval. 


i".  *  - 


•‘.V.V.V.V, 


Keystrokes 


Purpose 


Keystroke 
number  for 
first  stroke 
on  line 


000 

LBL  A 

002 

RCL  30  x  i  t 

put  value  from 

30 

into  test  register 

005 

RCL  02  ST0  01 

old  becomes 

new 

eo 

009 

RCL  04  ST0  03 

old  cti  becomes 

new 

“o 

013 

RCL  26  *  it  *  180  «  STO  02 

calculate  0i  in 

radians 

024 

RCL  27  *  it  t  180  »  STO  04 

calculate  cti  in 

radians 

035 

(RCL  02  -  RCL  01)  *  RCL  05  =  STO  06 

store  a  in  06 

048 

(RCL  04  -  RCL  03)  *  RCL  05  =  STO  07 

store  c  in  07 

061 

RCL  06  -  RCL  07  =  STO  08 

store  (a-c)  in 

08 

069 

RCL  06  +  RCL  07  =  STO  09 

store  (a+c)  in 

09 

077 

0.5  *  RCL  08  =  STO  10 

store  l/2(a-c) 

085 

0.5  *  RCL  09  =  STO  11 

store  l/2(a+c) 

093 

(RCL  08  *  RCL  05)  cos  *  RCL  10 
=  STO  12 

store  A  in  12 

107 

(RCL  09  *  RCL  05)  cos  *  RCL  11 
=  STO  13 

store  B  in  13 

121 

(RCL  08  *  RCL  05)  sin  *  RCL  10 
=  STO  14 

store  E  in  14 

135 

(RCL  09  *  RCL  05)  sin  *  RCL  11 
-  STO  15 

store  F  in  15 

149 

RCL  03  cos  *  RCL  01  cos 
=  STO  16 

store  P  in  16 

159 

RCL  03  cos  *  RCL  01  ^in  =  STO  17 

store  Q  in  17 

169 

RCL  03  sin  *  RCL  01  cos  =  STO  18 

store  R  in  18 

179 

RCL  03  sl£  *  RCL  01  sin  =  STO  19 

store  T  in  19 

RCL  08  x2  fk 


189 


calculate  I ( a— c) I 


Keystroke 
number  for 
first  stroke 


line 

Keystrokes 

Purpose 

193 

INV  x>t  E 

if  (a-c)  is  not  zero,  calculate  Ax 
and  Ay  using  the  following 
equations.  If  (a-c)  equals  zero  (is 
within  0.0001  of  zero)  calculate  Ax 
and  Ay  using  the  equations  at  label  E 

196 

RCL  16  *  (RCL  10  + 

-  RCL  13)  +  RCL  17 

-  RCL  18  *  (RCL  14 
*  (RCL  12  -  RCL  13 
=  STO  20 

RCL  11  -  RCL  12 
*  (RCL  14  +  RCL 

-  RCL  15)  -  RCL 

-  RCL  10  +  RCL 

15) 

19 

ID 

calculate  Ax  and  store  in  20 

254 

RCL  16  *  (RCL  14  - 
*  (RCL  12  -  RCL  13 

-  RCL  18  *  (RCL  12 

-  RCL  11)  +  RCL  19 
=  STO  21 

RCL  15)  +  RCL  17 
-  RCL  10  +  RCL  11) 

+  RCL  13  -  RCL  10 
*  (RCL  14  +  RCL  15) 

calculate  Ay  and  store  in  21 

312 

LBL  B 

314 

RCL  06  x2 

calculate  |a| 

318 

INV  x>t  D 

if  a  is  within  0.0001  of  zero,  Az  is 
calculated  using  the  equation  at  D; 
otherwise,  Az  is  calculated  using  the 
following  equation 

321 

((RCL  06  *  RCL  05  +  RCL  01)  ^in 
-  RCL  01  sin)  *  RCL  06  =  STO  22 

calculate  Az  and  store  in  22 

344 

LBL  C 

346 

RCL  20  SUM  23 

calculate  x  (=£Ax)  and  store  in  23 

346 

RCL  20  SUM  23 

calculate  x  (= £ Ax)  and  store  in  23 

350 

RCL  21  SUM  24 

calculate  y  (=£Ay)  and  store  in  24 

354 

RCL  22  SUM  25 

calculate  z  (=EAz)  and  store  in  2  5 

358 


RCL  20 


display  Ax 


Keystroke 
number  for 
first  stroke 

on  line  Keystrokes  Purpose 


361  LBL  E 

363  (RCL  06  *  RCL  05)  sin  x2  *  2  * 

RCL  06  =  STO  31 

380  (2  *  RCL  06  *  RCL  05)  sin  *  4  * 

RCL  06  =  STO  32 

398  RCL  05  *  2  *  STO  33 

405  RCL  16  -  RCL  19  =  STO  34 

413  RCL  17  +  RCL  18  =  STO  35 

421  RCL  34  *  RCL  31  +  RCL  35  *  RCL  32  + 

RCL  33  *  (RCL  17  -  RCL  18)  =  STO  20 

446  RCL  31  *  RCL  35  -  RCL  32  *  RCL  34  + 

RCL  33  *  (RCL  16  +  RCL  19)  =*  STO  21 

471  GTO  B 

473  LBL  D 

475  RCL  01  cos  *  RCL  05  =  STO  22 


for  | (a-c) |  B  0 

2 

sin  (aAs)  , 

store  - r - -  in  31 

2a 


store 


sin(2aAs ) 
4a 


in  32 


store  —  in  33 

store  (P-T)  in  34 
store  (Q+R)  in  35 

store  Ax  in  20 

store  Ay  in  21 

return  control  to  location  B 
routine  calculates  Az  for  |a 
Az  in  22 


484 


GTO  C 


return  control  to  location  C 


APPENDIX  B:  TEXAS  INSTRUMENTS  Tl-59  CALCULATOR  PROGRAM  FOR 
QUADRATIC-VARIATION  MODEL,  ALLOWING  As  TO  VARY. 


Notation  is  as  before,  plus 

PGM  xx  -  program  xx  in  Master  Library  Solid  State  Software  module. 

INV  SBR  -  inverse  subroutine  returns  control  from  a  program  (or  other 
subroutine)  to  the  calling  routine. 


This  program  uses  the  Master  Library  PGM  09,  Simpson's  Approximation 
(continuous)  to  evaluate  the  integral  that  arises.  Because  of  restrictions 
in  the  logic  of  PGM  09,  it  is  necessary  to  calculate  Ax  (or  Ay  or  Az,  but 
only  one)  down  the  entire  borehole,  and  then  change  the  argument  of  the  in¬ 
tegral  to  calculate  Ay  or  Az.  Running  time  for  the  program  varies  with  the 
number  of  intervals  (n)  chosen  for  Simpson's  approximation.  Typical  values 
for  n  and  running  time  t  are  given  below  (for  calculation  of  Ay  between  1900 
and  1905  m  depth  at  DYE-3). 

n  t(s) 

2  21 

4  27 

6  36 

10  50 

20  88 


We  find  n=*6  to  be  a  good  compromise  between  speed  and  accuracy.  Increasing  n 
from  6  to  20  gave  a  change  of  only  1. 56x10” 6%  in  Ay  at  1900-1905  m.  In  many 
cases  n=4  is  quite  accurate  (to  within  0.01%  or  better). 

Model  used: 

Integrate  from  0  to  As2  (from  yo  to  yi): 

y  -  ♦,  a 

2  2 

Asi  (yi-yp)  -  a32  (yp-y-i),  2  rAsi  (yi-yp)  +  As 2  (yp-y-i>| 

y  A8i  As2  ( As i  +  As£)  3  As^  As2  (Asi  +  A$2)  +  y0 


if  Asi  *  As2 


ryi  -  2  yo  +  y-i  i  2  .  ryi  -  y-i , 
l - tt~2 - J  s  +  l — — — i  s  +  yo* 

2  As  2 As 


15 


.n’.nV.V 


Set  calculator  in  radians  (press  rad). 

Enter  program  into  calculator. 

Store  8_i  (in  radians)  in  memories  11  and  12. 

Store  81  (in  degrees)  in  13. 

Store  <x_i  otQ  (in  radians)  in  14  and  15. 

Store  <*i  (in  degrees)  in  16. 

Store  Asi  and  As2  in  17  and  18,  respectively. 

Store  0  in  01  (lower  limit  of  integration  for  calculating  Ax,  Ay  or 
Az).  Store  n  (we  prefer  n=6)  in  28. 

Press  A. 

After  about  30  seconds,  calculator  displays  Ax  (for  Ay  or  Az,  depending  on 
the  argument  programmed  for  the  integral). 

Memory  29  contains  the  running  sun  of  Ax  (or  EAy  or  £Az). 

Program  automatically  reassigns  0o  ♦  6-i ,  01  +  +  °<-i  ,  «i  +  «0J 

so,  to  calculate  next  step,  simply  store  9i,  cq  (in  degrees)  for  new 
step  in  memories  13  and  16,  respectively,  and  press  A.  Note  that  if 
As i ,  As 2  change,  it  is  necessary  to  store  the  correct  new  values  in 
memories  17  and  18. 

Note  that  this  procedure  is  designed  to  work  down  (or  up)  the  hole  and 
will  not  apply  to  the  first  interval  at  the  top  (or  bottom)  of  the  hole. 

The  value  of  Ax  (or  Ay  or  Az)  can  be  calculated  using  the  linear-variation 
routine  and  can  be  stored  in  memory  29  so  that  memory  29  contains  the  true 
running  value  of  the  coordinate  being  calculated.  Alternately,  one  can 
reverse  direction  (if  calculating  downhole,  calculate  the  first  segment 
going  uphole)  and  store  the  negative  of  the  resulting  A  value  in  memory 
29.  Thereafter,  one  follows  the  instructions  at  the  top  of  this  page. 


Keystrokes 


Purpose 


Keystroke 
number  for 
first  stroke 
on  line 


000 

LBL  A 

002 

RCL  13  *  tt  t  180  =  ST0  13 

Change  0i  from  degrees  to  radians 

013 

RCL  16  *  tt  t  180  =  ST0  16 

Change  cq  from  degrees  to  radians 

024 

RCL  17  *  RCL  18  *  (RCL  17  + 
=  ST0  19 

RCL  18) 

040 

(RCL  17  *  (RCL  13  -  RCL  12) 
(RCL  12  -  RCL  ID)  *  RCL  19 

-  RCL 
=  STO 

18  * 
20 

069 

(RCL  17  *  (RCL  16  -  RCL  15) 
(RCL  15  -  RCL  14))  *  RCL  19 

-  RCL 
=  STO 

18  * 
22 

098 

(RCL  17  x2  *  (RCL  13  -  RCL  12)  + 
RCL  18  x2  *  (RCL  12  -  RCL  11))  * 
RCL  19  =  STO  21 

129 

(RCL  17  x2  *  (RCL  16  -  RCL  15)  + 
RCL  18  x2  *  (RCL  15  -  RCL  14)) 
t  RCL  19  =  STO  23 

160 

RCL  28  STO  05 

store  n  in  05  for  discrete 
integration 

164 

RCL  18  STO  02 

store  As z  in  02  for  discrete 
integration  (integrate  from  0  to  AsD 

168 

RCL  18  5-  RCL  05  -  STO  03 

calculate  h  (discrete  interval  width) 
and  store  in  03 

176 

PGM  09  D 

call  integration  program  and 
integrate 

179 

RCL  04  SUM  29 

update  running  sum  in  29  of  A  values 

183 

RCL  12  STO  11 

reassign  0g  +  0_i 

187 

RCL  13  STO  12 

reassign  -*■  @o 

191 

RCL  15  STO  14 

reassign  uq  *  a-i 

195 

RCL  16  STO  15 

reassign  oti  au 

199 


RCL  04 


display  calculated  &  value 


Keystrokes 


Purpose 


Keystroke 
number  for 
first  stroke 
on  line 


To  calculate  Ax 


202 

LBL  A' 

204 

((ST0  10  x2  *  RCL  22  +  RCL  10  * 

RCL  23  +  RCL  15)  cos  *  (RCL  10  x2  * 
RCL  20  +  RCL  10  *  RCL  21  + 

RCL  12)  sin) 

calculate 

Ax 

243 

INV  SBR 

return  control  to  main  program 

To  calculate  Ay  (Identical  to 

Ax  except  cos 

sin) 

202 

LBL  A' 

204 

((ST0  10  x2  *  RCL  22  +  RCL  10  * 

RCL  23  +  RCL  15)  sin  *  (RCL  10  x2 
*  RCL  20  +  RCL  10  *  RCL  21  + 

RCL  12)  sin) 

calculate 

Ay 

243 

INV  SBR 

To  calculate 

Az 

202 

LBL  A’ 

204 

((STO  10  x2  *  RCL  20  +  RCL  10  * 

RCL  21  +  RCL  12)  cos) 

calculate 

Az 

224 

INV  SBR 

18 
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